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Sum rules for the nonperturbative piece of correlators (specifically, the vector current correlator) are discussed. 
The sum rule subtracting the perturbative part is of the superconvergent type. Thus it is dominated by the bound 
states and the low energy production cross section. It leads to a determination of the gluon condensate (a B G 2 ). 
We find 

(q s G 2 ) ~ 0.048 ± 0.030 GeV 4 . 



1. SUM RULE 

The potential, or more generally the spectrum of 
a system of heavy quarks cannot be directly dis- 
cussed in terms of the OPE (operator product ex- 
pansion). However, one can use dispersion rela- 
tions to deduce a number of sum rules relating 
bound state properties to quantities obtainable via 
the OPE ("ITEP-type" sum rules). One can then 
use the estimates of nonperturbative contributions 
to bound states energies and wave functions to 
actually go beyond the traditional analysis. Al- 
though the sum rules, being global relations, can- 
not discriminate details one can check consistency 
and even obtain reasonable estimates on nonper- 
turbative quantities, specifically on the gluon con- 
densate. This last is the aim of the present note, 
where we will use a method generalizing that pro- 
posed by NovikovW. 

To do so we consider the correlator for the 
vector current of heavy quarks, 



IV =(P 9^ -PnPv) u (P ) 
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where = V'TmV' an d sum over omitted colour 
indices is understood. This will give information 
on triplet, I = states; information on states with 
other quantum numbers would be obtained with 
other correlators. The function II(£) satisfies a 
dispersion relation, 



n(t) = ^J d: 



P(s) 



t 



(2) 



where p(s) = Imil(s). Actually, this equation 
should have been written with one subtraction. 



We will not bother to do so as its contribution 
drops out for the quantities of interest for us here. 

Let us denote by n p t ., p p .t. to the cor- 
responding quantities calculated in perturbation 
theory, albeit to all orders, but nonperturbative 
effects are neglected in n p .t., p P .t.- (In actual cal- 
culations we cannot of course include all orders. 
We will sum the one-gluon exchange to all orders 
(which can be done explicitly in the nonrelativistic 
regime), and add one loop radiative corrections to 
this.) In particular, for example, the gluon con- 
densate contribution is not included in the "p.t." 
pieces. 

At large t, both spacelike and timelike, 
the OPE is applicable to il(t), and we have the 
well-known results! 2 ! 



n(0 = iW<> + ^ 



(3) 



and 



, v N C C F (a s G 2 > (l + v 2 )(l-v 2 ) 2 
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with v = (1 — 4m 2 /s) 2 the velocity of the quarks. 
Moreover, 
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N c = 3, C F =4/3. 
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If we then define IInp , Pnp as the results of sub- 
tracting the perturbative parts, 

n NP = II - ITp. t . ; p NP = p - pp. t ., 

it follows from the OPE, Eq. (3), that n NP de- 
creases at infinity like t~ 2 and hence satisfies a 
superconvergent dispersion relation. We thus have 
a first sum rule: 



/ 



dspNp(s) = 0. 



(5) 



In fact it would appear that one still has another 
sum rule because of the following argument. At 
large t, Hnp(£) behaves like (cf. Eq. (3)) 

while the contribution from the bound states to 
the dispersion relation (see below), 

IlNP;bound states (t) ~ ^2q,3 " ^) 

dominates over this. Therefore we have the extra 
relation, 



ds spjqp(s) = 0. 



(7) 



It turns out that (7) is actually equivalent to (5), 
up to radiative corrections. This is because the 
region where any of the integrals in (5), (7) are 
appreciably different from zero is for s ~ 4m 2 (1 + 
0(a 2 )), so (7) differs from (5) by terms of order 
a 2 , smaller than the radiative corrections which 
neither (5) nor (7) take into account. 

Let us return to the sum rule (5). The 
function p(s) consists of a continuum part, for s 
above threshold for open bottom production, and 
a sum of bound states. Both can be calculated 
theoretically provided that s is larger than a cer- 
tain critical s(v ), and n smaller or equal than a 
critical no- s(vo) and n are defined as the points 
where the perturbation theoretic contribution to 
p and the nonperturbative one are of equal mag- 
nitude, and form the limits of the regions where a 
full theoretical evaluation is possible. 

To be precise, for the continuum we use 
(4) so that above the critical s(vo), 

N C C F (a s G 2 ) (l + v 2 )(l-v 2 ) 2 



P c nT(s) 
s > s(v ), 
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(8a) 



and v is such that p^ , p t (s(v )) ~ /3 p ° t nt (s(u )); 
numerically, and for bb, vo ~ 0.2. For the bound 
states p is proportional to the square of the wave 
function at the origin, 

N 

p(s) = J ^\R n (0)\ 2 6(s-M n ). 

We may get Pp[l,(s) and p^p{s) by splitting the 
residue |i?„(0)| 2 into a Coulombic piece, 



I!?' 



Coul. 



(0)| 2 
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where the one loop corrections S^'Os may be 
found in ref. 3, and the (leading) nonperturbative 
correction are given by the Leutwyler-Voloshin 
analysis (refs. 4, 3). So we have 

|i?„(o)| 2 ^|^ ouL (o)l 2 + l^ oul -(o)l 2 C p ; 

the numbers 5^ P have been calculated by Leutwyler 
and Voloshin. For n = 1, 
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This is all we really need since, for bottomium, 
n = 1. Thus we have 



PnpM = 



3N c C 3 F Trm 3 (a s G 2 ) 



8a?m 4 



n=l 



n <no, 



(8b) 

the r] n known in terms of the 5^ P ■ 

The sum rule (5) can then be written 
schematically as 



/ Pnp + ^ Residue of p NP 



rs(v ) 



'threshold 



Pnp + ^ Residue of p^p > . 



n—UQ + l 



The left hand side is given in terms of (a s G 2 ) by 
Eqs. (8); the right hand side can be connected with 
experiment with the following argument. The sum 
over higher bound states, 

"E„=„ 0+ i Residue of p NP " 
may be identified as the difference between the 
sum over the experimental residues of the poles 
of the bound states, and what we would get by a 
Coulombic formula, for all n > hq + 1. Certainly, 
this Coulombic formula will not be valid for large n 
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because here the radiative corrections will become 
large; but, because the residues decrease like 1/n 3 
the contribution of these states will be negligible. 
We write this decomposition as 

(bound states with n > n ) =p^ P ,n>n ( s ) 

PCoulombic, n>n 

As for the continuum piece below s(vq) we may 
likewise interpret it as the difference between ex- 
periment and a perturbative evaluation, which we 
write as 

= PI^(S) P C X(S), * < 8{vo), 

and, because we are close to threshold, we have 
N c C F a s 1 



cont / \ _; 
F P-t- V 7 8 1 _ e -irC F a 3 /v 



(1 + ^cont^s) 



N c C F a s 



(1 + S cont a s ) 



and the value of the one loop radiative correction 
^cont^s may be found in ref. 5. 

Taking everything into account the sum 
rule (5) becomes, 

E :^K xp (o)l 2 + /bac k M 
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We have defined vq = ea s and the expression 
is valid up to corrections of relative order a s . 
The function /back(fo) is the contribution of the 
background which, when added to the resonances 
above threshold (included in the sum in the l.h.s. 

of (9)), give the experimental value of /^sLid Pnp- 
The function /back would be obtained by integrat- 
ing the cross sections for production of T + G and 
BB, where by G we mean a "glueball" decaying 
into 2tt, and B is any of the states B°,B ± ,B*. 
Because we may assume that the structure is pro- 
vided by the resonances, we can take /back given 
by phase space only. So we have 

/back(vo) = fl v /2 + h V l 

where the first term refers to the channel T + G , 
and the second to BB. We have in this expression 
neglected tog- 



2. NUMEROLOGY 

In principle the procedure would appear straight- 
forward. One would fit the resonance and bound 
state residues and / to the data, and then, af- 
ter substituting into (9), obtain a determination 
of (a s G 2 ). In practice, however, things do not 
work out so nicely. The quality of the experi- 
mental data does not allow any precise determina- 
tion of the constants /i^; any values in the range 
fi.2 ~ 0.03, 0.1 would do the job. Secondly, the 
effective dependence of (a s G 2 ) in Eq. (9) on ex- 
periment is proportional to aj 6 : so the result will 
depend very strongly on the value of a s we choose. 
This is particularly true because radiative correc- 
tions to the nonperturbative contribution to the 
bound states have not been calculated, so there is 
not even a "natural" renormalization point. 

These two difficulties may be partially 
overcome with the following tricks. First of all, 
since we are assuming that the n = 1 bound state 
is described with the bound state analysis as dis- 
cussed in ref. 3, we may fix the value of a s that 
produces such agreement. This means that we 
will take 0.35 < a s < 0.4. Secondly, we may alter 
the treatment of the continuum in the following 
manner. We split not from «o, but from vi, ar- 
bitrary provided only that v\ > v$. Thus, for 
s < s( Vl ), we use p^(s) = f%g(a) - p- t nt ( S ), 
and for s > s(vi) we take the theoretical expres- 
sion 



N C C F (a s G 2 ) {\+v 2 )(l-v 2 f 
128 s 2 v 5 



The sum rule is thus written as 

= 2C 3/ [c(3) _ lla »_i5^!) 



n=2 



48e 3 a 3 m 4 J ' 



eia s = Vi. 



Then we may profit from the fact that the sum 
rule should be valid for all values of v\ > vq to fix 
/i,2 requiring this independence, at least in the 
mean. That is to say, that when we increase v\ 
past a particle threshold from T(2) to T(6) the 
variation of the corresponding determinations of 



(a s G 2 ) around their average be minimum. The 
calculation may be further simplified replacing 

/backM -> 2f vf 75 . 

The results of the analysis are summarized in the 
following tables, where the column "Res" indicates 
at which resonance the cut in v\ occurs. We have 
taken two rather extreme values of /q. 



Res. 




(a s G 2 ) 


T(2) 


0.21 


0.014 


T(3) 


0.34 


0.034 


T(4) 


0.40 


0.048 


T(5) 


0.43 


0.039 


T(6) 


0.46 


0.046 



For a s = 0.35, / = 0.04 



Res. 


Vl 


(a s G 2 ) 


T(2) 


0.21 


0.037 


T(3) 


0.34 


0.057 


T(4) 


0.40 


0.067 


T(5) 


0.43 


0.048 


T(6) 


0.46 


0.052 


For 


a s = 0.40, f 


= 0.09 



This derivation shows very clearly the kind 
of errors one encounters. To the variations that 
may be called "statistical" , apparent in the differ- 
ent values found in the tables above 

0.014 < (a s G 2 ) < 0.067 

we have to add "systematic" ones, e.g., the influ- 
ence of the not calculated radiative corrections, 
easily of some 30%: not to mention our includ- 
ing the Coulombic wave functions at the origin for 
large values of n, or the lack of definition of the 



expression "perturbation theory to all orders" be- 
cause of renormalon ambiguities. Given all these 
uncertainties, which do even make it dubious that 
one can really define with precision the conden- 
sate in terms of experimental observables, it is 
not surprising that one cannot pin down the gluon 
condensate with more accuracy than an estimate, 
taking into account above figures, of 

(a s G 2 ) ~ 0.048 ± 0.03 GeV 4 . 

To get this average we have taken into account all 
determinations in the tables above, excluding the 
lowest (T(2)) and highest, T(6). This is slightly 
larger than old averages, and slightly lower than 
more recent ones I 6 ' which tended to give, respec- 
tively, (a s G 2 ) ~ 0.042, (a s G 2 ) ~ 0.065 GeV 4 . 
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